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Abstract Alfven waves play three related roles in the impulsive phase of a solar 
flare: they transport energy from a generator region to an acceleration region; 
they map the cross-field potential (associated with the driven energy release) 
f->i • from the generator region onto the acceleration region; and within the acceler- 

ation region they damp by setting up a parallel electric field that accelerates 
electrons and transfers the wave energy to them. The Alfven waves may also be 

^O I regarded as setting up new closed current loops, with field-aligned currents that 

close across field lines at boundaries. A model is developed for large-amplitude 
Alfven waves that shows how Alfven waves play these roles in solar fiares. A 

(VS ' picket-fence structure for the current fiow is incorporated into the model to 

^ I account for the "number problem" and the energy of the accelerated electrons. 

Mh. Keywords: solar fiares; Alfven waves; electron acceleration 

o 

C/2 ' 1. Introduction 

^. 

There are long-standing, unsolved problems in the physics of solar flares. A large 

fraction of the magnetic energy released in a flare appears in e = 10 — 20 keV 
^ • electrons that produce hard X-rays and type III solar radio bursts. One problem 

r^ ! is that there is no satisfactory model for the acceleration of these electrons, 

Qv ' which, in the older literature, was referred to as "flrst phase" acceleration (Wild, 

Smerd, and Weiss, 1903) and as "bulk energization" of electrons. Runaway ac- 
celeration (Holman, 1985) due to a parallel electric field [_E||] seems the only 
f~^ [ viable acceleration mechanism, but there is no accepted model for how i?|| ^ 

ff^ ' is set up and maintained in an acceleration region. There is also a long-standing 

"number problem" (Hoyng, Brown, and Van Beck, 1976; MacKinnon and Brown. 
1989; Bian, Kontar, and Brown, 2010) that can be expressed in various ways; for 
, , , example, the rate, N > 10^^ s~^ , of precipitation of accelerated electrons inferred 

rS ' from hard X-ray observations implies a total of 10'^^ accelerated electrons in a 

jrt ! flare of duration 10^ seconds, and this exceeds the number of electrons stored 
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Figure 1. A cartoon showing a model due to Fletcher and Hudson (2008) for energy release 
in a flare in which energy is transported from an energy-release site (effectively at the top 
of the loop), via a propagating twist to an acceleration region in the chromosphere. In the 
Aifven model in this article, the flux loop is replaced by a vertical magnetic field, with the 
energy-release site at the top. 



in the entire flaring flux loop, estimated to be 10'^^ (Ein.slic and Henoux, 1995). 
Another problem is that the power released in a magnetic explosion (Melrose, 
2012a) may be expressed as the product [/$] of a current and a potential, and if 
one assumes a current of order / = 10^"'^ A, consistent with vector magnetogram 
data, then to account for the energy released one requires a very large potential, 
$ = 10^° V, for which there is no direct evidence. The latter two problems are 
related, in that there is a single large unexplained factor [M] with eN = MI 
and 6$ — Me, suggesting M of order 10^. To account for the factor M, Holman 
(1985) argued that there are multiple current paths, with the current [/] flowing 
up and down M times. However, no mechanism has been proposed to account for 
this seemingly bizarre current pattern. The resupply of electrons is attributed to 
a return current (Brown and Bingham, 1984; Spicer and Sudan, 1984; van den 
Oord, 1990; Litvinenko and Somov. 1991; Emslic and Henoux. 1995). Arguments 
related to the number problem and the return current suggest an acceleration 
region near (or in) the chromosphere, where there is an adequate supply of 
electrons. 

A model developed by Fletcher and Hudson (2008) addresses the number 
problem. The model is shown in Figure 1. Energy release in the corona due to 
magnetic reconnection generates large-scale Aifven wave pulses, which propagate 
within a loop to the chromosphere where they damp and accelerate electrons 
via a turbulent cascade. The model includes some of the essential features in 
an acceptable flare model. Notably the energy-release site and the accelera- 
tion/dissipation region are remote from each other, with transport of energy 
between them via Aifven waves. 

The model envisaged here (Melrose, 2012b) for energy release involves a 
generator (or energy-release) region, where reconnection allows conversion of 
inflowing magnetic energy into outflow in the form of plasma kinetic energy and 
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an Alfven flux (c/. Fletcher and Hudson, 2008). We do not discuss the details 
of the generator region here, but several comments are appropriate. The energy 
release needs to be driven, and one possible mechanical driver is a pressure 
gradient, associated with collection of dense plasma near the top of the flaring 
loop (Haerendel, 2001, 2009; 2012, Fletcher and Hudson, 2008). In this case 
one needs to consider the energy and momentum associated with this dense 
plasma, and how it is transferred to Alfven waves. Our view (Melrose, 2012a, 
2012b) is that the driver is the Maxwell stress, and that this becomes available as 
reconnection allows the magnetic figuration to change to a less stressed state. The 
energy and momentum in the plasma then play only a minor role. We identify 
the power [/$] released to redirection of pre-existing current across field lines in 
the generator region, with the electromotive force [<J>] attributed to the rate of 
change of the stored magnetic flux. This cross-field current [/] is assumed to be 
redirected so that it flows along field lines through an acceleration region, and 
closes across field lines in the photosphere, and then flows back to the generator 
region. 

A schematic diagram of the new model is shown in Figure 2. The magnetic 
fleld is assumed to be vertical, with the top of the figure in the corona, and the 
bottom at the photosphere. The redirection of current at the coronal generator 
is analogous to the redirection of current that occurs in the so-called current 
wedge in a magnetospheric substorm (McPherran, Russcl. and Aubry, 197.3; 
Paschmann, Haaland, and Treumann, 2002). The physics underlying this effect 
is related to a conducting boundary (Simon, 1955), which is the plasma wall 
in a laboratory context, the ionosphere in the magnetospheric context, and the 
photosphere in the context of a flare. Specifically, when a current is driven across 
field lines in the body of the plasma, it closes by flowing along the field lines 
to a region where the conductivity allows it to close across field lines. Following 
Holman (1985), it is assumed that this redirection occurs M times in a flare, 
so that the potential imposed at the generator region is $/Af for each such 
redirected current path. 

In this article an analytic model for the large-scale Alfven waves is developed 
and used to describe various roles that Alfven waves play in such a model. The 
features that are new, in the solar context, are the transport of potential along 
fleld lines, and the setting up of new closed current loops. The electric field in 
the wave is expressed in terms of a wave potential, which also determines the 
field-aligned current (FAC) density in the wave. The waves map the potential, 
imposed at the generator region, onto the acceleration region, where the waves 
are assumed to damp, converting the cross-field potential at the top of the 
acceleration region into a field-aligned potential within the acceleration region. 
Figure 3 illustrates the mapping of the potential and the role of the acceleration 
region. With the magnetic field in the vertical direction, the equipotentials shown 
by dashed lines are vertical above the acceleration region. The potential set up 
in the generator region (above the top of the diagram) maps along field lines 
to the acceleration region, where a field-aligned potential is established. The 
field-aligned electric field [Eu] in the region extracts energy from the waves and 
transfers it to electrons. 

Assuming that current is conserved in an Alfven wave (V • J = 0), an Alfven 
wave can transfer a FAC from one field line to another, and can transfer a 
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Figure 2. A schematic diagram of the model adopted in this article, as explained in the text. 



cross-field current from one location to another, but it cannot introduce any 
intrinsically new current. This implies that an Alfvtsn wave that is localized to a 
flux tube can have no net FAC: the FAC in the wave must change sign across the 
wavefront, and include equal up and down currents. The up and down FACs in 
the wave close across field lines in the generator region and in the photosphere, 
forming a closed current loop, and transferring the cross field current from one 
region to the other. 

The general model for Alfvcn waves is developed in Section 2, and some 
relevant boundary conditions are discussed in Section 3. The model is applied to 
cylindrical and planar geometries in Section 4. Some aspects of the extension 
of the model into the acceleration region arc discussed in Section 5, and a 
general discussion of the results is presented in Section 6. The conclusions are 
summarized in Section 7. 



2. Alfven Waves 

The model for Alfven wave adopted here is based on one developed in connection 
with magnetospheric substorms (Vogt and Hacrcndcl, 1998; Paschniann, Haa- 
land, and Trcumann, 2002). Conventional treatments of Alfven waves are based 
on MHD or the low-frequency limit of cold-plasma theory, and assume that the 
waves are of small amplitude. The approach adopted here applies to Alfven waves 
even when neither of these theories is valid, and there is no restriction on the 
wave amplitude. Our approach is based on two well-known results from orbit 
theory: the electric and polarization drifts (Northrop, 1963). It allows one to 
include the plasma response to an electric field (and its time derivative) without 
assuming that its amplitude is small. 
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2.1. Electric and Polarization Drifts 

Any electric field, including both the inductive electric field due to the changing 
magnetic field, and the electric field in an Alfvcn wave, can be separated into 
components perpendicular [E^] and parallel [£'||z] to the background magnetic 
field [B = Bz] which is assumed uniform along the z axis. Subject to relatively 
weak conditions (on the slowly time- varying electric field), one has 



E — El 



S||Z, 



E^ = u X B, 



(1) 



with u = Ej_ X B/S^ the electric drift velocity. In MHD, u is interpreted 
as the fluid velocity, and Equation (1) is attributed to Ohm's law for infinite 
conductivity, which also implies E^\ = 0. The interpretation of u as the electric 
drift velocity is different from the usual MHD interpretation: here E implies u 
and i^ii is arbitrary, whereas in MHD u implies E with E\\ = 0. The divergence 
of the electric field, in the form Equation (1), gives the charge density. 



— = -Bz-uj 



dE, 



dz 



a; = V X u. 



(2) 



where uj is the vorticity in a fiuid interpretation. 

A temporally changing (perpendicular) electric field corresponds to a dis- 
placement current e{)d¥\±/dt, and this causes a polarization drift, which is 
along d'E±/dt at a velocity proportional to the mass to charge ratio. After 
summing over the contributions of all species of charged particles, this implies a 
polarization current density 



3l± 



-I 



displi 



'displ 






(3) 



The sum of the two currents in Equation (3) is (? jv^ times Jdispi, where wq = 
?;a/(1 + '^\I<?Y^'^ is the MHD speed and wa is the conventional Alfven speed. 

Equation (3) may also be derived using cold-plasma theory, as shown in 
Appendix A. 

2.2. Wave Equation 

Including the polarization current in Ampere's equation, and combining it with 
Faraday's equation, leads to a wave equation for the electric field. The parallel 
component of the cold-plasma response is given by Equation (42), and including 
it in the parallel component of the wave equation gives 



V^E, 



zV^^i 



Vp 



= V 



Vp 1 a2E_L . 




Id-EA 


c2 at2 


--I) 


V_L-E_ 


dz 





(4) 



(5) 
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For Eh ^ the term involving the charge density couples the perpendicular and 
parallel components in Equation (4). 

The wave equation Equation (4) simplifies for E\\ = and V_l x E^ = 0, 
with the latter condition implying V^(V_l • E^) = V^^Ej^. Then Equation (4) 
becomes the usual wave equation 






Ei = 0. (6) 



2.3. Wave Potential 

A general solution of Equation (6) is of the form 

Ej_(a;, y, z, t) = E+(a;, y, z + vqI) + E2(a;, y,z~ VQt). (7) 

Assuming the z axis is vertical, the + solution is propagating down, and is 
referred to as a direct wave, and the — solution is propagating up, and is referred 
to as a reflected wave. The condition Vj^ x E^ = must be satisfied for Equation 
(7) to be valid, and this implies that Ej^ may be written 

E± = -Va^^^. (8) 

iterpretcd 
advection equations: 



The quantity $_,_ is interpreted as the wave potential. The two solutions satisfy 



|Ei . ±„„|Ei. (9) 

Using Equation (9), the magnetic and electric fields in the wave are related by 

woB^ = ±z X E^. (10) 

The electric drift velocity may be interpreted as the fluid velocity [u^ = Ej_ x 
B/B^] in the wave. It satisfies the Walen relation in the form 

^^A (U) 

vo B 

where Equation (10) is used. 

The energy density in the waves satisfies the well-know equipartition relation 
for Alfven waves, 

i£o|E±|^ + i^|u±p = ^, (12) 

where ?/ is the mass density. 

2.4. Parallel Current 

It is convenient (Song and Lysak, 2000) to define an effective current density [J'] 
that is the sum of the actual current density and the displacement current, 

aE 

J' = J + eo^, V-J' = 0, (13) 
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where the latter fohows from MaxwcU's equations. The perpendicular current 
associated with the waves is J'^ = J'j^ + J'^ with 

J'± = ^^-±^^ (14) 

where Equation (9) is used, and where Ra = Mo^o is the Alfven impedance. One 
has 



which implies 



J,f = T^V,.Ei^±J-Via>i. (16) 



Assuming E^^ ~ 0, Equation (16) determines J|| = J!.. 



3. Boundary Conditions 

The Alfven wave model shows how energy is transported, potential is transferred 
and new current loops are set up. In applying this to an idealized flare model, 
one needs to impose boundary conditions in three regions: the generator region, 
the acceleration region, and the photosphere. Prior to a flare, the acceleration 
region does not exist, and its turning on is identified with the triggering of the 
energy release in a flare. 

3.1. Reflection Coefficients 

An Alfven wave incident on a boundary between between two plasmas can be 
partly reflected and partly transmitted. The "waves" here have effectively infinite 
wavelength, in the sense that the boundary regions may be regarded as of zero 
thickness compared with the wavelength. 

In the magnetospheric application, the lower boundary is the ionosphere, 
which is partially ionized, implying nonzero Pedersen and Hall conductivities. 
Let the reflection coefficient be fref- The relation between the reflected and 
incident waves is (Scholcr, 1970) 



Ej_ = r,.ef E;};, J|| =-r,.ofJ|p r^cf = - p , p , (17) 



Ra + Rp 



where 1/Rp is the height-integrated Pedersen conductivity. The reflection co- 
efficient is zero in the case of perfect absorption, and this corresponds to the 
impedance matching condition Rp = Ra- 

In the solar context, the chromosphere and photosphere are partially ionized, 
and one could model the effect of the photosphere using Equation (17) with 
Ra ^ Rp- A more relevant effect is the increase in density from above to 
below the photosphere/chromosphere, implying a decrease in Alfven speed. For 
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long wavelength waves, the reflection coefficient at a boundary between two 
non-dissipative plasmas, labeled 1 and 2, is 

^ Rai - Ra2 ,^n-. 

Rai + Ra2 

with Rai,2 = y'OVAia- 

The simplest approximation to the photospheric boundary condition corre- 
sponds to line-tying, which is the limit of infinite inertia and zero Alfven speed, 
Ra2 ^' in Equation (18). Thus the line- tying boundary condition corresponds 
to Ticf — > —1. Line-tying implies E]] = — EJ[, u~ = — u"*", B^ = Bj], JJ = J^, 
Jr = J|t. Hence there is no electric field or fluid velocity at the boundary, 
but there is a net current and associated magnetic field. A stress imposed in 
the generator region is transferred to the photosphere by the waves, where the 
resulting J x B force is balanced by the assumed infinite inertia. This stress is 
transferred back to the generator region by the waves. 

One also needs a boundary condition at the generator region. Suppose this 
is Tref = —1, which corresponds to a constant- voltage generator (Paschmann, 
Haaland, and Treumann. 2002). At such a boundary, the electric fields in the 
direct and reflected waves cancel [E^ + EJ = 0] and the parallel currents add 
[J^ = Jr] with this net current modifying the cross-field current in the generator 

region. The opposite assumption [r^ef = 1] corresponds to a constant-current 
generator. This boundary condition corresponds to 'E']_ = — E]] and jJ' + J7 = 0, 
so that the cross-field current in the generator is unchanged, and the electric field 
and flow velocity are modified. Intermediate cases between these two limits have 
been discussed in the literature on substorms (Lysak, 1985; Vogt, Haerendcl, 
and Glassmeier, 1999). 

Prior to a flare, the stresses are in balance, field lines are equipotentials and 
there is no energy transport. Any stress imposed on the generator region is 
prevented from causing the plasma to move (dynamically) by line-tying in the 
photosphere. 

3.2. Boundary Condition on the Flux Tube 

Alfven waves transfer current across field lines, allowing redistribution of field 
aligned current (FAC), but they cannot generate intrinsically new FACs. The 
redistribution of FAC may be interpreted as a current loop associated with the 
Alfven waves, with up and down FACs closing across field lines in both the 
generator region and the conducting boundary. The net FAC current associated 
with the Alfven waves, found by integrating J^ + J7 over the cross-sectional 
area of the fiux tube within which the waves are confined, must be zero. Using 
Equation (8) and Equation (16), this leads to a restriction on the form of $^ -t- 
$^. Assuming that this restriction applies separately to ^^, the condition is 



n-V$=^L = 0, (19) 



where S is the surface of the fiux tube and n is the normal to it. For example, 
in a cylindrical model for a flux tube of radius i?. Equation (19) requires that 
d$±/dr = at r = i?. 
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4. Cylindrical and Planar Models 

Cylindrical and planar models for the system of Alfvcn waves arc developed in 
this section. For simplicity, impedance matching is assumed, so that there is no 
reflected wave. 

4.1. Cylindrical Model for the Potential 

A polynomial model for the potential that satisfies Equation (19) is 

d$1 



- ± 



dr 



Ar^iR-r)", (20) 



with A, a, b constants. It follows from Equation (16) with Equation (20) that one 
has 



J|l =-^Ar''-\R-r)''-^[{a + l)R-{a + b+l)r]. (21) 



1 

rI' 

To avoid singularities at r = and r ~ R one needs a > 1 and b > 1, respectively. 
The sign of the FAC reverses at r = ro = (a + l)R/{a + 6+1), such that the 
direct and return currents flow at r < vq and r > 7'o. This dividing radius [rg] 
decreases with increasing b. One may integrate Equation (20) to find an explicit 
expression for $+(r). For a > 1 and b >1, the potential difference $+(r) — $+(0) 
does not change sign in the region < r < R. 

The Poynting flux at a given r is proportional to (d^^/dr)^, which goes to 
zero for r ^f R. The total energy flux is determined by integrating the Poynting 
vector over the flux tube. The direct and reverse currents are found by integrating 
Jii over the ranges < r < ro and r^ < r < R, and are equal in magnitude. Both 
the total energy flux and the direct and return currents increase with increasing 
&, corresponding to decreasing rg. 

4.2. Planar Models for the Potential 

Let the spatial variation perpendicular to the magnetic fleld be along the x 
axis, with no dependence on y. It is convenient to assume that the Alfvcn wave 
propagates between boundaries at x = xi and x ~ X2. 
A polynomial model that satisfies Equation (19) is 

'^ = A[x- (xi + X2)l2r[{x - Xi){x2 - x)]^ (22) 

ax 

where A, c, and d are constants. This model is similar to the cylindrical model, 
with 

J\\ = -f^{c{x ~ x{){x2 -x)- 2d[x - {xi + X2)l2f], (23) 

which has one sign near the center of the flux tube, at a; = {xi + X2)/2, and the 
opposite sign nearer the edges at x = xi and x ^ X2- The sign of J\\ reverses at 
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x±, with 



x± 



Xi + X2 



± 



2c X2 — Xi 



2c + d 



(24) 



The energy transport and the direction and return current in this model are 
quahtatively similar to those in the cylindrical model Equation (20). In particu- 
lar, both increase with increasing d in Equation (22), which plays a similar role 
to b in Equation (20). 

The planar model facilitates incorporating many pairs of direct and return 
currents. For example, consider the sinusoidal potential, which satisfies the con- 
dition Equation (19), 



d$+ 
dec 



= 5 cos 



X - {xi +a;2)/2 

X2 - Xi 



Mtt 



with b a constant and M an integer. This model gives 



bMir 



Ra{x2 -Xi) 



X - {Xl +X2)/'2. 
X2 - Xl 



Mtt 



(25) 



(26) 



In this case it is straightforward to integrate Equation (25) to find the potential 



*+ A+ b{x2~xi) . 
<I>t - $n = ^ sm 



X - {Xi+ X2)/2 
X2 - Xl 



Mtt 



Ra{x2 -Xifj\\ 

M^TT^ 



(27) 



where ^q is a constant. For M = 1, the direction of the FAC for xi < x < 
(xi + X2)/2 is opposite to that for (xi + X2)/2 < x < X2, and one has a single 
pair of direct and return FACs. The inclusion of the multiplicity factor [M] 
implies M neighboring pairs of up and down FACs. 



5. Inclusion of Dissipation 

The foregoing discussion applies to undamped Alfven waves propagating between 
the generator region and the acceleration region. Once the waves enter the ac- 
celeration region it is assumed that they are damped through some collisionless 
dissipation process. A self-consistency argument allows one to relate the E\\ ^ 
involved in the acceleration to the spatial decay of the waves. 

5.1. Onset of Effective Dissipation 

The onset of a flare is attributed to the turning on of effective dissipation in 
the acceleration region, due to the turning on of E\\ ^ 0. The dissipation 
must be anomalous: this point has been recognized by many authors (Alfvcn 
and Carlcivist, 1967; Colgate, 1978; Holman, 1985; Melrose and McClymont, 
1987). Although various models for the anomalous dissipation have been explored 
(Raadu, 1989; Borovsky. 1993; Bian. Kontar. and Brown, 2010; Hacrcndcl. 2012), 
no consensus has emerged on the detailed microphysics involved in the creation of 
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i?ll ^ 0. An approach that avoids these details is based on the idea that the rate of 
anomalous dissipation is able to adjust to meet global requirements (Haercndel, 
1994; Melrose, 2012a). All of the proposed forms of anomalous dissipation imply 
highly localized and transient regions of dissipation, and effective dissipation 
requires a statistically large number of such localized transient regions coupled 
together. On a macroscopic scale, the effective rate of dissipation is determined 
by the statistical distribution of these localized, transient regions, and is in- 
sensitive to the details of the microphysics (Melrose, 2012a). In the present 
context, this implies that the rate of dissipation due to acceleration by -Eii in 
the acceleration region can adjust to the rate of dissipation required by the 
rate of magnetic energy release in the energy-release region. We assume that 
this adjustment is relatively slow, over many Alfven propagation times, in a 
precursor phase. We identify the onset of a flare as a threshold being reached 
where the coupling between the localized dissipation regions becomes effective, 
Uke a percolation threshold in a network system (Dorogovtsev and Mendez, 
2002). In such a model, the macroscopic dissipation turns on and couples to 
the generator region in at most a few Alfven propagation times. The power 
dissipated in the acceleration region can then be modeled as Regl^, where ReS 
is an effective resistance, which adjusts to meet the required rate of dissipation. 
The development of En ^ implies that the frozen-in condition does not 
apply in the acceleration region, allowing the magnetic field lines to slip through 
the plasma. This was referred to as a "fracture" of the magnetic field by Haercn- 
del (1994; 2012). Whereas prior to the flare, the line-tying condition precludes 
plasma motion in response to the stress imposed in the energy-release region, 
once i?|| j^ develops in the acceleration region, the plasma above this region 
can move relative to the plasma below it. This motion corresponds to the fluid 
motion [u^] in the Alfven waves, and is essential in allowing Alfvenic trans- 
port of energy. The assumption that the potential surfaces all close within the 
acceleration region is an oversimplification, that ignores the requirement that a 
return current be set up between the acceleration region and the photosphere, or 
more specifically, the region of cross- field current closure. As in the application to 
auroral acceleration (Marklund, 2009), this requires some S-shaped, as well as U- 
shaped potentials. Once effective dissipation turns on, the rates of energy release 
and dissipation adjust to balance each other within a few Alfven propagation 
times. 

5.2. Impedance Matching 

The foregoing arguments suggest that the global requirement on the dissipation 
and electron acceleration is that Alfvenic flux released in the energy-generation 
region be completely absorbed in the acceleration region. This condition corre- 
sponds to impedance matching. 

Suppose the dissipation is described by an effective (anomalous) conductivity, 
aes, which is averaged over the distribution of localized dissipation regions. Let 
1/RoS be the effective conductivity integrated over height in the acceleration 
region. The condition for zero reflection of an Alfven wave incident on the 
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Figure 3. The lines of constant potential are along vertical field lines above the acceleration 
region, and are shown to close as semicircles within the acceleration region. The potential 
gradient changes from horizontal to vertical with decreasing height. 



acceleration follows by analogy with Equation (17): 

Rcs — Ra 



rrct 



RcS + Ra 



0. 



(28) 



It follows that all of the incoming energy is dissipated completely when the 
effective resistance of the acceleration region adjusts such that it is equal to the 
Alfvenic impedance. This impedance matching condition [i?off = ^a] leads to a 
model for the dissipation that is independent of the details of the microphysics 
involved in the anomalous dissipation. 

5.3. Dissipation and the Parallel Electric Field 

Outside the acceleration region, one can assume E^ = in the waves, so that 
the field lines are equipotentials, and the frozen-in condition applies. Within 
the acceleration region E'^ ^ implies that neither of these conditions apply. 
The plasma slips through the magnetic field, such that all fields associated with 
the wave, E]^, $][, etc., decrease with decreasing z. In particular, this implies 
that $]| has a dependence on z that is not included in the dependence through 
z ± vot. This leads to a component of the potential electric field along z: 



e;^ 



OZ Vq otf 



(29) 



This electric field is assumed to accelerate particles, transferring energy from the 
waves to the particles, leading to the damping of the waves. 

The conventional interpretation of how En ^^ is set up was suggested by 
Gurnett in 1972 (Alfven, 1977), and is illustrated here in Figure 3. The idea is 
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that the cquipotcntial surfaces, which are along field lines above the acceleration 
region, must close across field lines within the acceleration region. Then E'J' is 
attributed to the normal to cquipotcntial surfaces having a component along 
field lines. A physical argument is that the cquipotcntial surfaces must close 
above a perfectly conducting boundary, and that the fluid velocity associated 
with the waves must go to zero above the line-tied boundary; that is, both E_l 
and uj^ must go to zero above the photosphere. The following analytic model 
illustrates how this leads to i?|| 7^ 0. 

In a cylindrical model, above the acceleration region the field lines are equipo- 
tentials, so that $* depends only on r and z±vot. Within the acceleration region, 
it is assumed that the cquipotcntial surfaces close across the field lines. Suppose 
the acceleration region is at 2 < 24. The shape of an equipotential surface at 
z < Zt can be described as a function of r,z. Let the potential surfaces be 
described by the equation $(r, z) — $0 = f{i'iz) = 0, where different surfaces 
depend (implicitly) on a constant r^, that is the value of r at z > z^. The ratio of 
the parallel and perpendicular components is then determined by the function 

The lower boundary of the acceleration region is the surface /(r, z) — Q for 
rt = R, where rt = R is the boundary of the cylinder at z = zj. A particle 
that passes through the acceleration region along a given field line at radius r is 
accelerated through a potential difference 

A$+ = |$+(i?) -$+(r)|, (31) 

where the right hand side is to be evaluated at z > Z4. 

A specific model for the surfaces is to assume that they are elliptical. One has 
vertical equipotential surfaces [r = rt < R, for z > zt] and elliptical surfaces, 

/(r, z) = r^+ a" {zt - zf - r^ = 0, (32) 

for z < zt^ with a the ratio of the senfi-axes of the ellipse. The case a = 1 
corresponds to semi-circular surfaces illustrated in Figure 3. The lower boundary 
of the acceleration region is the surface r^ -\- a'^{zt — z)^ — R^ =0. The ratio 
Equation (30) becomes 

E+ ■.E+^a^{z-Zt):r. (33) 

The boundary condition Equation (19) implies that both Ef and E\_ vanish on 
(and below) the lower boundary of the acceleration region. 

While plausible, this model is heuristic. The assumptions made is writing 
down the wave solution Equation (7) and introducing the wave potential Equa- 
tion (8) arc not valid for E]\ ^ 0. 
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6. Discussion 

The three features of the Alfven wave model emphasized here are energy trans- 
port, the mapping of potential differences, and the setting up of intrinsically new 
current loops involving FACs. In this section, the cylindrical and planar models 
described in Section 4 are used to illustrate these features. 

6.1. Alfvenic Energy Transport 

Energy transport by Alfven waves is well understood: energy propagates at 
vq « va along the direction of the background magnetic field. Such energy 
propagation has been invoked in several different flare models. For example, 
Piddington (1974) developed a model in which the Alfven waves transport the 
energy into the corona from below the photosphere during the flare. In the flare 
model due to Fletcher and Hudson (2008), and in the model developed here, 
energy already stored in the corona is converted into an Aflvenic energy flux 
in a generator region, is transported downwards, and is transferred to energetic 
electrons in an acceleration region. 

In the Alfven wave model, the energy transport is closely linked to the FAC 
[J||] in the waves. It is of interest to compare the present model for energy trans- 
port with an earlier model (Melrose, 20121)) in which, prior to the flare, the FAC 
profile in a cylindrical flux tube is described by a function ji(^), with ^ = r/R. 
This profile changes, to J2{S,), after passage of an Alfvenic front launched by the 
onset of the flare. It was found that power transported increases with increasing 
concentration of J2(^) towards ^ = 0; this increase continues as the central con- 
centration of J2(0 increases to arbitrarily large values, compensated by a return 
current flowing at larger radii. In the present model, the pre-flare FAC is excluded 
by the assumption that the field lines are straight. This may be interpreted as 
assuming that the guiding magnetic field is arbitrarily strong, and that both the 
pre-flare FAC and the FAC in the Alfven waves may be treated as (independent) 
perturbations. The Alfven wave model then has the same qualitative and semi- 
quantitative features as in the earlier model. Specifically, in the cylindrical case, 
the power transported by the Alfven waves increases as Jn becomes increasingly 
concentrated towards the axis of the cylinder, where it enhances the pre-flare 
FAC, with J|| in the opposite direction to the pre-flare FAC at larger radii. In 
particular, for the model Equation (20) both the concentration of the FAC and 
the energy flux increase with increasing power-law index b. The planar model 
Equation (23) has similar properties, with the power-law index d the counterpart 
of b in the cylindrical model. 

6.2. Transport of Potential 

A new feature of the Alfven wave model in the solar context is the transport 
of potential. The Alfven waves are launched in the generator region, where the 
power released is due to a cross- fleld potential and a cross-field current: the 
cross-field potential in the wave matches that at the boundary of the generator 
region, and the FAC in the wave arises from redirection of the cross-fleld current 
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generator 



A$ = ^tot/M 



aromiloui coiductivity' 



A$ 



photosphere 



Figure 4. The model envisaged here has the generator region at the top, with a total potential 
<E> and cross-field current /, belowr which is an energy transport region to an acceleration region 
with anomalous conductivity. The cross-field current is redirected along field lines to close in 
a conducting boundary, which is assumed to be the photosphere. Following Holman (1985), 
the current [/] is assumed to flow up and down multiple [M] times, forming a picket-fence 
structure. The potential across the acceleration region is A<I> = <I>/M. 



at the boundary of the generator region. Which of these two boundary conditions 
dominates determines whether the generator can be approximated as constant- 
vohage or constant-current, which in turn depends on the ratio of the internal 
and external impedances. The Alfven waves transport the potential along field 
lines to the acceleration region. Reflection of the Alfven waves from both regions 
modifies the potential, providing a feedback that regulates the rate of dissipation 
in the acceleration region and the rate of energy release in the generator region. 

A realistic theory requires a physical model to determine the form of the wave 
potential. Here, simple analytical models are chosen to satisfy the boundary 
condition Equation (19) at the edge of the flux tube within which the Alfven 
waves are confined. The magnitude of the potential is not constrained by the 
theory, but it is implausible that the total potential available [^tot — 10^ — 
10^*^ V] can be transported by a single Alfven wave. There is no evidence that 
^tot appears across the acceleration region; the evidence is that the potential 
that does appear across the acceleration region is A$ = $tot/-^, with M of 
order 10^. 

The explanation for the multiplicity [M] proposed by Holman (1985) can be 
incorporated into the model. Holman's suggestion is that there are M pairs of 
up and down current channels, leading to the "picket- fence" model illustrated in 
Figure 4. One requires that the potential across the acceleration region [A$ = 
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^tot/M] be of order 10^ V for each channel, each of which has up and down 
currents [A/] of order /. It is unreahstic to use cyhndrical geometry to describe 
muhiplc pairs of up and down currents, which would occur in alternate concentric 
rings. In a planar geometry, the pairs of up and down current are in alternate 
current sheets. A more realistic geometry requires an arrangement for the M 
current pairs, in either cylindrical (r, </>) or cartesian (x, y) geometry, combined 
with a model for each of the current pairs. 

A physical explanation is needed for the multiplicity [M] of current pairs. 
One possibility is to consider restrictions imposed by the microphysics on the 
potential in the Alfven waves. If the potential is restricted to $+ < 10"* V, 
then M =10^ current channels must develop in order to provide the required 
total rate of energy loss. Moreover, these current channels must develop si- 
multaneously; if they develop sequentially, the implied timescale of 10^ Alfven 
propagation times greatly exceeds the timescale of a flare. These problems are 
not discussed further here. 

6.3. Dissipation and Acceleration 

How Alfven waves damp and accelerate particles effectively is an unsolved prob- 
lem. Some form of anomalous dissipation is required. As has long been recog- 
nized, a requirement for a current-driven instability is a very high current density, 
which implies filamentation of the current into many high-current-density chan- 
nels (Holman, 1985; Tsuneta, 1985; Melrose and McClymont, 1987; Melrose, 
1990; Emslio and Henoux, 1995; Tsuneta, 1995). Let the current threshold for 
instability be Jcrit- Suppose each filament is cylindrical, and of radius A, such 
that the current in each filament is If = ttA^ Jdit- One requires a large number 
[iVf = I /If] of such channels. For a specific model, Hacrcndcl (2012) esti- 
mated Jcrit = 4 X 10'^ A m^^ and A = 1 meter. An unavoidable implication is 
that anomalous dissipation requires that the dissipation region is highly struc- 
tured, involving very small scales perpendicular to the field lines. The following 
argument based on the Alfven wave model leads to a similar conclusion. 

Suppose one balances the incoming Poynting vector in the Alfven waves with 
the outgoing kinetic energy flux in accelerated electrons with a number density 
Tig. This gives 

1 (<i^Hr)\'_^^ ( ^\''\e^^^rr^ (34) 



/iQi'o \ dr J \m 

with A$+(r) given by Equation (31), and where [2eA$+(r)/m]^/^ is the speed 
of the accelerated electrons. Using Equation (31), one may integrate Equation 
(34) to find 

eA$+(r)=2m.o2(J)', i = ^, (35) 

and use Equation (16) to find 

3 / r \ 2 
|J||| = en^vo- [jj . (36) 
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With eA$+(r) of order 10'^ x tuvq under coronal conditions, it follows from 
Equation (35) that the model requires perpendicular structure in $+(r) on a 
scale of a few tens of skin depths [c/wp] which is of the same order of magnitude 
as the A estimated by Hacrcndel (2012) by a different argument. It follows from 
Equation (36) that there is a high current density, sufficient to trigger anomalous 
conductivity, on the scale A. 

The result. Equation (36), demonstrates an inconsistency in the Alfven model, 
as developed here. The Alfven wave model applies on a macro scale, plausibly 
describing the energy and potential transport to the top of the acceleration 
region. Effective dissipation involves micro-scale structures, and these must be 
taken into account within the dissipation region. The model developed here for 
the dissipation region can at best be regarded as applying to properties averaged 
over a statistically large distribution of such micro-scale structures. 

6.4. Return Current 

The long-standing "number problem" can be partly resolved by requiring a re- 
turn current between denser regions of the solar atmosphere and the acceleration 
region. Existing models for the return current have invoked both electrostatic 
and inductive effects (Brown and Bingham, 1984; Spicer and Sudan, 1984; van 
(Icn Oord, 1990). The Alfven wave model provides a new way of modeling the 
return current, in terms of Alfven waves setting up new current loops, that 
close across field lines in the acceleration region and in denser regions of the 
solar atmosphere. Such a model has an important qualitative difference from 
these earlier models, in which the ions were neglected (van don Oord, 1990). 
Neglecting the ions leads to neglecting the polarization current, so that the 
inductive effects propagate effectively at the speed of light. It is essential to 
include the polarization current, and then Equation (3) implies that inductive 
effects are transported along field lines by Alfven waves at the MHD speed. An 
Alfven wave model for an inductively driven return current is needed to show 
how the electrons are continuously resupplicd to the acceleration region, as is 
necessary for both "first phase" electron acceleration and for auroral electron 
acceleration (Marklund, 2009). 



7. Conclusion 

The main point of this article is that large-amplitude Alfven waves play an 
important role in solar flare physics. These are not "waves" in the conventional 
sense, with a well-defined frequency and wave vector, but are specific solutions 
of the Alfven wave equation Equation (6) involving a cross-field wave potential 
[^j_] and field aligned currents, J,, , determined by the wave potential. Energy 
transport is along field lines at the MHD speed [vq ~ wa] as for any (torsional) 
Alfven wave. The wave potential provides the mapping of the potential imposed 
in the generator region onto the acceleration region. Reflected Alfven waves 
transport the reaction (of the acceleration region or the photosphere) back to 
the generator region, providing a coupling between them. The Alfven waves set 
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up new closed current loops, with oppositely directed FACs closing across field 
lines at the two ends (in the generator region and the photosphere). 

A long-standing problem with models involving particle acceleration by Alfven 
waves is the mechanism that allows the Alfven waves to damp and to transfer 
their energy to particles. The interpretation of the wave amplitude in terms 
of a cross-field potential provides a way of interpreting this damping without 
identifying the microphysics involved. Within the acceleration region it is as- 
sumed that the wave amplitude, and hence the cross-field potential, decreases 
downwards. There is then a nonzero gradient of the potential along the field lines 
implying a nonzero i?|| . This -En accelerates the particles. In an idealized model, 
the Alfvenic energy flux at the top of the acceleration region is converted into 
a kinetic-energy flux in accelerated electrons at the bottom of the acceleration 
region, where the wave amplitude (potential) is zero. The rate of dissipation can 
be determined by assuming that the anomalous dissipation processes adjust to 
maximize the power dissipated. This leads to the impedance matching condition 
for the effective resistance, Rcs = Ra, of the acceleration region. All anomalous 
dissipation processes involve microphysics on tiny space scales, of order one meter 
in the present context (Haerendel, 2012), and the relation between processes on 
micro and macro scales is treated heuristically here. This relation needs to be 
explained in a more detailed model. 

The Alfven wave model indicates how long-standing, unsolved problems re- 
lated to bulk energization of electrons may be resolved. It is assumed that the 
acceleration results from E\\ ^ in an acceleration region, along field lines 
above a hard X-ray source in the chromosphere. One problem is the seeming 
inconsistency between the large cross-field potential, of order 10"'^° Volts, and the 
energy, eA$ of order lO"* eV, of the accelerated electrons. One way of resolving 
this (Holman, 1985) is to assume that the pre-flare current flows back and forth 
between the generator and dissipation regions many (M = 10^) times. Electrons 
are accelerated downward in a sequence of M jumps in potential, one jump 
for each instance of / flows up through the acceleration region. The additional 
insight that the Alfven wave model provides is that M current paths may be 
regarded as M closed current loops set up by Alfven waves. The number problem 
is resolved by identifying the total rate [N] that electrons precipitate as M times 
the rate [I/e] for each of these loops. However, there is no obvious explanation 
within the model for the value of the multiplicity [M] . A speculation is that the 
microphysics in the acceleration region constrains the value of eA$. 

There are unsolved problems related to the acceleration region (Haerendel, 
2012). In particular, the self-consistency argument for E^^ ^ 0, based on the 
damping of the waves due to energy transfer to particles by E\\ , is heuristic, and 
needs to be complemented by a speciflc mechanism that accounts for E\\ ^ 0. All 
such mechanisms appear to involve very small scales perpendicular to the field 
lines, estimated to be of order one meter by Haerendel (2012). An outstanding 
problem is how this microphysics relates to the macroscopic model discussed in 
this article. 

The Alfven wave model developed here is highly simplified, and is applied 
specifically only to transport of energy between the generator and acceleration 
region. However, the physics involved is likely to be widely relevant to energy 
transport and energy release in the solar atmosphere. 
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Appendix 

A. Response of a Plasma at Low Frequencies 

The response of a cold plasma may be described by the dielectric tensor Kijiuj) 
(Stix, 19G2): 

/ S -iD \ 

Dduj) = eoIu.j{uj)Ej{uj), K,j{i^)=[iD S . (37) 

\ OP/ 

At sufficiently low frequencies, when dissipation is neglected, one has 



S'«l + ^, DwO, P=l--^. (38) 



The electric induction [D] includes the response of the plasma through the 
polarization [P]: 

ap 

D = £oE + P, pi„d = -V-P, Ji„d = 7^, (39) 

at 

where pi^d and Ji„d are the induced charge and current densities. The perpen- 
dicular component of the response is 

Pi = — eoE±. (40) 

The temporal derivative of Equation (40) gives 

2 nxp 

which reproduces Equation (3). The parallel term in the low- frequency, cold- 
plasma limit gives 

The components Equation (41) and Equation (42) of the current density are 
included in Maxwell's equations. For the perpendicular components, one obtains 

where the right hand side includes source terms. The result Equation (43) re- 
produces Equation (6) when the source terms are neglected. For the parallel 
component, one obtains 



1^^^_V72 



1 9pcxt ^-'"extll .... 
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